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Abstract 

We consider a special type of triple pendulum with two pendula attached to end mass 
of another one. Although we consider this system in the absence of the gravity a quick 
analysis of of Poincare cross sections shows that it is not integrable. We give an ana- 
lytic proof of this fact analysing properties the of differential Galois group of variational 
i i equation along certain particular solutions of the system. 

> 1 Introduction 

in 

In this paper we study dynamics of a special type of triple pendulum which we called 
a flail. It consists of three pendula. The first one is attached to a fixed point, and to its 
end mass the other two pendula are joined. The lengths and the masses of pendula are 
hf hf h> an d ni\, m 2i m 3f respectively, see Fig. [T| We consider this system in the absence 
^vq of gravity. Similar problem of dynamics of a simple triple pendulum in the absence of 

t-H gravity field was analysed in [TSJ. 

The configuration space of our system is a torus T 3 = S 1 x S 1 x S 1 with local coordi- 
nates (9\, 8 2 , 8 3 ) mod 2n. We chose the point of suspension of the first pendulum as the 
origin, and angles are measured from vertical line as it is shown in Fig. [T] The Lagrange 



X 



function has the following form 

L = l 



* (l 2 l} iel + l 2 2 m 2 6 2 2 + Z|m 3 0|) + hl 2 m 2 cos(0 1 - e 2 )^6 2 

+ hi 3 m 3 cos(e 1 -e 3 )0 1 e 3r 



where \i = m.\ + m 2 + m 3 . Let us remark that the dynamics of this system can be in- 
terpreted as a geodesic motion on torus T 3 . The system has S 1 symmetry. This is why 
the Lagrange function depends on differences of angles only. Hence, it is reasonable to 
introduce new variables defined by 

71 = 02-01, 72 = 03-02, 73 = 01 (1-2) 
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Figure 1 : Geometry of the flail pendulum 



(1.3) 



and corresponding momenta 

n\ = (l%m 2 + l 3 m 3 )ji + l 3 m 3 j 2 + [Zf w 2 + '3^3 + hhm 2 cos 71 + l\l 3 m 3 cos(7i + 7 2 )]7 3 , 
tt 2 = l 3 m 3 ji + l 3 m 3 j 2 + hm 3 [h + h cos(jt + j 2 )]j 3 , 

n 3 = [1%m 2 + l 3 m 3 + hl 2 m 2 cos 71 + hl 3 m 3 cos(7i + 7 2 )]7i + hm 3 [l 3 + h cos^ + 7 2 )]72 
+ \i\m 2 + \\m 3 + + 2/!/ 2 m 2 cos 71 + H\l 3 m 3 cos(7i + 7 2 )]7 3 . 

Now, 73 is a cyclic variable, and 7T3 is a first integral. Thus, after this reduction, the system 
has two degrees of freedom and its dynamics is determined by Hamiltonian 

H = [ll]i{llm 3 {Kx - n 2 f + \\m 2 T^ 2 ) + /^m^Oi - c) 2 
- l x (l x l\m\n\ cos7f + / 3 m 3 cos(7 1 + 72)(2/ 2 ?« 2 7r 2 (c - n x ) 
+ l 1 l 3 m 3 (n 1 - 7r 2 ) 2 cos(7i +72)) - H 2 l 3 m 2 m 3 (n x - 7r 2 )cos7i(/3(7r 1 - c) 
+ l\Tt 2 cos(7! + 72)))] / [2llljl 2 m 2 m 3 B], 

where 

B :— }i — m 2 cos 2 71 — m 3 cos 2 (7i + 72), 

and c : = 713 is a parameter. 

The aim of this paper is an analysis of integrability of this reduced system. To get 
an idea about the complexity of dynamics we made several Poincare cross sections. In 
Figures|2j|3j and|4]we show such cross sections. We take 71 = and tz\ > as the cross 
section plane. The phase portraits are presented in (72, y 2 ) plane. The cross sections are 
made for the following values of parameters: 

7t 3 = c = 1, m 1 = l, m z = 3, m 3 =2, l x = 1, l 2 =2, l 3 = 3. (1.4) 

Figures are ordered according to increasing energy of the system. They show that, for 
chosen values of parameters, the system is not integrable. The main problem considered 
in this paper is to prove that for a wide range of the parameters the system is in fact 
non-integrable. Our main result we formulate in the following theorem. 



Theorem 1.1. If the reduced flail system governed by Hamiltonian 1 1.3 1 has parameters satisfying 

• m\ 7^ 0, c 7^ 0, l 2 7^ l 3 , and m 2 l 2 = m 3 l 3 ; or 

• m 2 = m 3 , l 2 = l 3 , and c = 0, 

then it is nonintegrable in the class of meromorphic functions of coordinates and momenta. 
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Figure 2: The Poincare sections for £ = 0.01 (on the left) and £ = 0.01005 (on the right) 
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Figure 3: The Poincare sections for £ = 0.0102 (on the left) and £ =0.011 (on the right) 
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Figure 4: The Poincare sections for £ = 0.012 (on the left) and £ = 0.013 (on the right) 
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The plan of the paper is the following. Section [2] contains the non-integrability proof 
and is divided into two subsections. In Section|3]some remarks are given. Explicit forms 
of some complicated expressions are given in Appendix [A] For convenience of readers 
in Appendix [B] basic information about linear second order differential equation with 
rational coefficients and Kovacic algorithm are given. 

2 Non-integrability proof 

In this section we prove Theorem |1 . 1 | using the Morales-Ramis theory, see lffl [TTl . 

We formulate the main theorem of this theory in the following simple settings. Let 
us consider a complex Hamiltonian system. We assume that the underlying phase space 
is C 2 ", which is considered as a linear symplectic space equipped with the canonical 
symplectic form 

n 

U) = Yj dqi A dp ; , 

i=l 

where x = {cj\, ...,q n ,pi,...,p n ) are the canonical coordinates. Let H : C 2n — >■ C be a 
holomorphic Hamiltonian, and 

^x = v H (x), v H {x)=l 2n V x H, xeC 2n , feC, (2.1) 

the associated Hamilton equations. Here I^m denotes 2n x 2n symplectic unit matrix. 

Let t i — y <p(t) 6 C 2n be a non-equilibrium solution of ( |2.1) . The maximal analytic 
continuation of <p(t) defines a Riemann surface T with f as a local coordinate. Variational 
equations along q>(t) have the form 

^£ = A(f)£, A(t) = d ^( <P (t)). (2.2) 

We can attach to this equations the differential Galois group Q. Roughly speaking, Q is 
an algebraic subgroup of GL(2n, C), which preserves polynomial relations between so- 
lutions of | |2.2) . As a linear algebraic group, among other things, it is a union of a finite 
number of disjoint connected components. One of them, containing the identity, is called 
the identity component of Q, and is denoted by Q°. For a precise definition of the differ- 
ential Galois group and differential Galois theory see e.g., l5ll2l lT0llT4l . It appears that 
the integrability of the considered system manifests itself in properties of the differential 
Galois group of the variational equations. For Hamiltonian system and integrability in 
the Liouville sense this relation is particularly elegant. For these systems Q is an alge- 
braic subgroup of Sp(2n, C). In nineties of XX century Morales-Ruiz and Ramis showed 
that the integrability in the Liouville sense imposes a very restrictive condition on Q° 

nuinnni. 

Theorem 2.1 (Morales-Ruiz and Ramis). Assume that a Hamiltonian system is meromorphi- 
cally integrable in the Liouville sense in a neigbourhood of a phase curve I. Then the identity 
component of the differential Galois group ofNVEs associated with I is Abelian. 

In order to apply this theorem we need an effective method which allows to deter- 
mine properties of the differential Galois group of linear equations. In the investigated 
system variational equations split into two subsystems of linear equations. Each of these 
subsystems can be transformed into equivalent second order equation with rational co- 
efficients. For such equation there exists the Kovacic algorithm 0, see also Appendix|Bj 
which allows to determine the differential Galois group. We apply this algorithm in our 
considerations. 
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2.1 Case m-i ^ m.3 and I2 7^ I3 

At first we make the following non-canonical transformation 



"?r 






72 


= A 








2/1 


. 7T 2. 




J/2. 



A = 



1 
-2 








1 1 - 
1 



(2.3) 



After this transformation equations of motion take the form 



where 



J 2 „:=BIB\ B = A 



H(x) := H(Ax). 



(2.4) 
(2.5) 



The explicit form of Hamiltonian in these variable is following 



H 



\ 2 m 2 m 3 (ci 3 + hyz - h{y\ +yi)) + Hf '2^23/2 + mihyi - hyi) 



- h (m 3 cos(x 1 - x 2 ) {2l\m 2 y 2 (c/ 3 + Z2J/2 - '3(1/1 + J/2)) + hm 3 (hyi ~ hyi f cos(x! - x 2 , 
+ 2l 2 m 2 m 3 {-l 3 yi + l 2 y 2 ) (-l 2 y 2 + l 3 (-c + y x + y 2 ) + hy 2 cos{x l - x 2 ))cosx 2 
+ l\l 2 m 2 y 2 cosx 2 \ / [2l 1 l 2 l 3 m 2 m 3 (}i — m 3 cos(xi — x 2 ) — ffi2COsx 2 )]. 

(2.6) 

An explicit form of vector field v(x) := S 2n X7 x H is given in Appendix |a| see equa- 
tion ( |A.1| ). These equations are still Hamiltonian but written in non-canonical variables. 
Let us assume that the following conditions 



m 2 l 2 = m 3 l 3 , l 2 ^ l 3 , 
are fulfilled. Then, system <|2.4| has an invariant manifold 



(2.7) 



Af = {{xi,x 2 ,yi,y 2 ) G C 4 | x 1 = y 1 = 0, y 2 = -cl 3 /{l 2 - l 3 )}. 

In fact, restricting the right hand sides of system | |2.4| to J\f , and putting m 3 = m 2 l 2 /l 3 , 
we obtain 

±1 =0, x 2 = - - n rT , y 1 =0, 1/2 = 0. (2.8) 



l 2 m 2 (l 2 -l 3 )' 

Hence, solving the above equations we obtain our particular solution C 3 1 1— > <p{t). 

Let X := [Xi,X 2/ Yi,Y 2 ] T denote variations of [x\,x 2 ,yi,y 2 ] T . The variational equa- 
tions along cp(t) take the form 



df 



X = AX, 



A = 



a n a 13 a u 

a 21 a 23 «24 

«3i a 33 fl 34 

#41 #43 #44 



(2.9) 



The explicit forms of entries are given in equation \A.2\ in Appendix [Aj 

We note that equations on X\, Y\ and Y 2 do not depend on X2, i.e., variable X2 decou- 
ples from remaining variables. Moreover one can observe that equations for Y\ and for 
Y 2 are dependent, see appropriate coefficients a,-,-. This fact is connected with the exis- 
tence of a first integral of variational equations. The expansion of energy integral along 
the particular solution gives the first integral of variational equations 



h := gradH(<p) -X 



l 2 l 3 m 2 (l 3 - l 2 ) 



-hYl + (h + h)Y2\. 



(2.10) 
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At the level h = we have 



h + h~ 



(2.11) 



Then, substituting this into the sub-system of variational equations for X%, Y\ and Y 2 and 
eliminating Y\ we obtain the following second order linear equation for X\ 

X x + aXi + = 0, (2.12) 

with coefficients 

a = [4c(Z 2 + Z 3 )(2Z 2 / 3 + h(h + h) cos x 2 )(l 1 l 3 m 1 + k(l 2 + h)m 2 + 2l 2 l 3 m 2 cosx 2 )smx 2 ] 
/[h(h ~ h){tt\l\m 2 + l\(l 2 + J 3 )(Z 3 mi + (l 2 + l 3 )m 2 ) +4hhh(l 2 + l 3 )m 2 cos x 2 )(2l 3 m 1 
+ (h + h)m 2 - (h + h)m 2 cos(2x 2 ))], 

b = [2c 2 (2/ 2 / 3 + h(l 2 + l 3 ) cos x 2 )(4hl 2 l 3 (l 2 + h)m 2 + (4Z 2 Z 2 m 2 + l\{l 2 + l 3 ){l 3 m x 
+ (h + h)m 2 )) cos x 2 )]/[hll{l 2 - hfm 2 {ti\\\m 2 + l\{l 2 + fe)(Z 3 f»i + (Z 2 + Z 3 )m 2 ) 
+ 4/i/ 2 / 3 (/ 2 + l 3 )m 2 cosx 2 )(2l 3 m l + (/ 2 + l 3 )m 2 - (l 2 + h)m 2 cos(2* 2 ))]. 

In order to transform this equation into an equation with rational coefficients we use the 
following transformation 

t I — > z = cosx 2 (f). 



(2.13) 



This change of variable together with transformations of derivatives 



d 
df 



^dz' 
■ sin x 2 x 2 ■■ 

- cos x 2 x\ 



d2_ 



(z) 



,d 2 



dz 2 
csin x 2 
Z 2 m 2 (Z 2 - l 3 )' 
2 

c z z 

"/2m 2 (/ 2 -/ 3 )2' 



'dz' 



(z) 2 



c 2 (l 



l\m\{\ 2 -hr 



convert equation (2.12 1 into 

X"- 



(z)X' + ^(z)X = 0, X = X lf 



(2.14) 



where the prime denotes the derivative with respect to z. The explicit forms of coefficients 
p(z) and q{z) are the following 



-(z + za) 



2az 



4(7 



z 2 -l (z 2 -l)a-aj A + a(4z + a- 
(crz + 2)(4cr + z(4 + a(a + olj))) 



ccj)' 



z^ (z 2 - 1) [d"(z 2 - 1) - aj) (4 + 0"(4z + c- + 0:7)] ' 
where we introduced notation 



h' 



h 

w 
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m 2 



o-:=cc + p. 



These parameters take non-negative values, and we assume that they do not vanish si- 
multaneously. 

Now, we make the classical Tschirnhaus change of dependent variable, 



X = zt>exp 

which transforms 1 2.14) into its reduced form 

r(z)w, r{z) 



w 



(2.15) 



(2.16) 
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with coefficient 



r(z) 



3_ 

16 L( Z 

ct 3 ^ 2 (5z 2 



-1)- 
«5 2 



(z-S) 2 (z + S) 2 (z- P y 



■ 4) + Acr 2 z (4 + ^ 2 - z 2 ) + Act (20 + <5 2 - 5z 2 ) + 64z 



32t7 2 (z 2 -l)(z 2 -^ 2 )(z-p) 



where 



*2 ._ 



1 



1 



(2.17) 



ocy 1 
^' p:= -^-4- 
Let us notice that coefficient r(z) depends only on two parameters 5 and cr. We assume 
that both are positive and S > 1 . 

If 6 7^ 1, and c<!> ^ 2, then equation < 2.16) has six pairwise different regular singular 
points 

Zl,2 = ±1, Z3,4 = ±<$, Z5 = P, Z 6 = 00. 

The respective differences of exponents at these points are following 



Ai = A 2 



1 

2' 



A 3 = A 4 = 0, A 5 = 2, 



1. 



Now, we can prove the following result. 
Lemma 2.2. For 

3^=1, and 5cr ^ 1, 
the differential Galois group of equation < 2,16) is SL(2, C). 



(2.18) 



(2.19) 



Proof. We prove this lemma by a contradiction. Thus, let us assume that Q 7^ SL(2, C). 
Then, according to Lemma B.l Q is either a finite subgroup of SL(2, C), or a subgroup 
of dihedral group or a triangular group of SL(2, C). 

At first we show that the first two possibilities do not occur. As A3 4 = 0, local so- 
lutions in a neighbourhood of z* = Z3, and z* = Z4, contain a logarithmic term. More 
precisely two linearly independent solutions W\ and of [ 2.16) in a neighbourhood of 
z* 6 {z3,Z4}, have the following forms 



W\{z) = (z - z*) p /(z), 102(2) = ^i(z) ln(z — z*) + (z — z*) p /i(z), 



(2.20) 



where /(z) and /i(z) are holomorphic at Z*, and /(z*) 7^ 0. The local monodromy 
matrix corresponding to a continuation of solutions along a small loop encircling z* 
counterclock-wise gives rise to a triangular monodromy matrix 



-27ri 
-1 



for details, see |9J. A subgroup of SL(2, C) generated by M* is not finite, and thus also 
the differential Galois group Q of ( 2.16) is not finite. As is not diagonalisable, Q is not 
a subgroup of dihedral group. Thus, Q is a subgroup of the triangular group of SL(2, C). 

In order to check if it is really the case we apply the Kovacic algorithm, for details see 
AppendixjB] According to it, if Q is a subgroup of the triangular group of SL(2, C), the 
first case of this algorithm occurs, i.e., the equation admits an exponential solution. 

To check this possibility we determine sets of exponents 



E 1 :={(l±A 1 )/2} / for i = l 6. 



Thus we have 
Ei = E 2 = 



4'4 



£3 = £3 



E 5 



2' 2 



E 6 = {0,1}. (2.21) 
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Next, according to the algorithm we look for elements e = [e\, e^, £3, £4, £5, £(>) or Cartesian 
product £ = nf=i Ef such that 



d(e) :=e 6 -J^ £,• 6 N , 



(2.22) 



/-l 



where No is the set of non-negative integers. In our case, there exists only one £ £ £ with 
this property, namely 

1111 1 1 S 



4' 4' 2' 2' 



(2.23) 



for which d(e) =0. Thus, according to the algorithm, we look for a polynomial P ^ of 
degree d(e), such that it is a solution of the following equation 



P" + 2coP' + (a/ + co 2 - r)P = 0, 



(2.24) 



where 



to 



5 

E 

!=1 



In the considered case, we have P = 1, so equation 1 2.24) gives equality 



(crz + 2) (cr 2 S 2 +4)z+4cr 
which cannot be fulfilled for a > 0. This end the proof. 



0, 



□ 



Now, let us analyse the cases excluded in Lemma 2.2 Case a = is equivalent to 
ol = I1/I3 = and j6 = = 0. So, it is equivalent to the case of zero length of the 
first pendulum. In effect, the system consists of two simple independent pendula and is 
obviously integrable. 

In the case of confluence of singular point Z3 and Z5 we prove the following. 



Lemma 2.3. Let us assume that 



6=->l. 



(2.25) 



Then the differential Galois group of equation 1 2.16) is SL(2,C). 

Proof. Then coefficient r in equation ( |2.16| simplifies to the following form 

3 3 1 1 

r ~ ~ 16(z-l) 2 ~ 16(z + l) 2 ~ 4(z-z 3 ) 2 ~ 2p(z-z 4 ) + 
3-p 3+p 1 



(2.26) 



16(p + l)(z + l) 16(p-l)(z-l) 2p(p2-l)(z-z 3 )' 



where Z3 = — Z4 = J := 7.1 v. Now, if <5 > 1, then equation | |2.16) has five singularities 
zi t 2 — ±1, Z3 and Z4, and Z5 = 00. 

One can easily check that the difference of exponents at z* = Z3 vanishes. Thus the 
local solutions around this point have the form 1 2.20) . One can repeat reasoning from 
the proof of Lemma 2.2 In order to check whether the first case of the Kovacic algorithm 
can appear at first we calculate the sets of exponents £, for i — 1, ... ,5 



£1 = E? = 



1 3 
4' 4 



E3 = 



E 4 = {0,1}, E 5 = {0,1}. 



(2.27) 



One can easily check that for each e = [e\,e%, £3, £4, £5) 6 nf=i E; we have 



d(e)=e 5 -J^e j £ N . 

i=l 



(2.28) 
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Thus, if 5 > 1, according to Kovacic algorithm, the differential Galois group of 1 2, 16) is 
SL(2,C). 

For 8 — 1, i.e., for a — 2, coefficient r simplifies to 

r < z > = "iRW (Z29) 

Now, equation 1 2.16| has only three regular singularities, i.e., it is the Riemann P equa- 
tion. Using the well known Kimura theorem [6J, one can check that its differential Galois 
group is SL(2, C). □ 

Remark 2.4 An equivalent form of condition aS = 2 isfollozving 

l\{h + h){hm 2 + him + m 2 )) - tt\%m 2 = 0. 

For nti = and Zj ^ 0, using the Kovacic algorithm, one can prove only that if the dif- 
ferential Galois group of variational equations along the mentioned particular solution is 
not SL(2, C), then it is a finite subgroup of SL(2, C). Calculations in first and second case 
of this algorithm are simple but in the third case becomes really complicated. Poincare 
sections given in Fig.[5]suggest that also in this case system is not integrable. 




m 1 *tfg 



el- 



Figure 5: Poincare sections for m\ = 0, m 2 ~ 3, W23 = 2, l\ = 1, l 2 = 2, Z3 = 3, 7T3 = c = 1, and for 
£ = 0.0095 (on the left) and E = 0.011 (on the right). Cross-section plane 71 = 



The proof of the first part of Theorem 1.1 is an immediate consequence of Lemma 2.2 
I2.3l and Theorem 12. II 

2.2 Case I3 = I2 and 111.3 = m 2 

If masses and lengths of two pendula attached to the first one are equal, the one can 
expect that behaviour of such system is regular. But Poincare sections suggest that also 
in this case the system is not integrable. 

In order to prove that the system is not integrable we apply again the Morales-Ramis 
theory. For Z3 = I2 and = m.2, the reduced Hamiltonian \2.6\ and the correspond- 
ing equations of motion $2A\ simplify considerably. In this case, for c = 0, the system 
possesses the following invariant manifold 

TV = {(x 1 ,x 2 ,y 1 ,y 2 ) e C 4 | x x = \j\ = 0}, 

and its restriction to this manifold takes the form 

*i = 0, *2 = -^/ 3/1=0, y 2 = 0. (2.30) 
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Figure 6: Poincare sections for m\ — 1, mi — = 2, l\ = 2, I2 = h = 1, /T3 
(on the left) and E = 0.0363 (on the right). Cross-section plane 71 = 



c = \ and E = 0.0035 



These equations determine our particular solution cp(t). Variational equations j2.2\ along 
this particular solution are also simpler. Namely, entries fl, 4 of matrix A vanish for 1 < 
i < 4. 

Now, equations for Xi and Y\ form a closed subsystem of variational equations. From 
them we obtain a second order equation for X\ . It takes form 1 2.12) with coefficients 



4i/2('2 + h cosx 2 ){li(m\ + 2m2) + 21inticosxi) sinx2 
l\(2l\mi + \\{m\ + 2mi) + 4ZiZ2 m 2COSX2)(wii + mi — W2 cos (2*2))' 
2y\{h + h cos xi)(4l 1 l 2 mi + (2l\mi + l\ (m x + 2m 2 )) cosx 2 ) 
l\l\mi^2\\mi + l\(m.\ + 2wi2) + Alilinii cos xi) (m\ +m 2 — m 2 cos(2x2)) 



Changing independent variable according to ( 2,13) we obtain equation of the form | |2.14[ | 
with coefficients 

-8z 2 a + 16z 4 a- 4a(2 + 7) + 6z 3 (2 + a 2 (2 + 7)) - z(6 + 7) (2 + a 2 (2 + 7)) 
V ~ (z 2 - l)(-2 + 2z 2 - 7) (2 + 4za + a 2 (2 + 7)) ' 

2(1 + za)(4a + z(2 + a 2 (2 + 7))) 



<7 



(z 2 - l)a(-2 + 2z 2 - 7) (2 + 4za + a 2 (2 + 7)) ' 
The coefficient r in reduced equation < 2.16) has the following expansion 
3 3 113 



+ 



+ 



16(z-l) 2 16(z + l) 2 4(z-<f) 2 4(z + <5) 2 4(z-z 5 ) 
; + a(15 + S 2 + a(6 + 5 2 {2 + a(-l + S 2 )))) 

16a(-l + S 2 )(1 + a(2 + a£ 2 ))(z - 1) 
8 + 3a(-5 + 2a) + a(-l + a (2 + a))^ 2 - a 3 <5 4 -2<S + a(-5 + £(<5 + a(-2 + aJ(-l + <5 2 )))) 



16a(-l + <5 2 )(1 + a(-2 + oc5 2 )){z + I) 
25 + oc{-5 + 5{5 + K{2 + ct5{-l+5 2 )))) 
4a<5(-l + 0i8) 2 {-\ + 5 2 ){z + S) 



4<x5(1+k5) 2 (-1+S 2 )(z-5) 
a (5 + 3a 4 <S 4 + oc 6 5 6 + a 2 (-16 + 75 2 )) 
-1 + a 2 <S 2 ) 2 (l + a 4 <5 4 + 2a 2 (-2 + 6 2 )){z - z 5 ) 

(2.31) 



where Z5 = — (1 + a 2 <5 2 )/ (2a) and S 2 = 1 + 7/2. Equation \2.14\ with coefficient r given 
in | |2.31) has singularities Z\ r2 = ±1, Z3 4 = ±<5 and Z5 = — (1 + a 2 5 2 ) / (2a), that are all 
poles of the second order and infinity has degree 3 for all positive a and 5. 
At first we consider the generic case. 
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Lemma 2.5. If 5 7^ 1/ol and 5^1, then the differential Galois group of equation [ 2,16) with 
r(z) given by < |23T) is SL(2, C). 

Proof. We apply the Kovacic algorithm. The sets of exponents E; in this case are 

~,|}, £ 3 = £ 4 ={^}, £ 5 ={~'|}' £ 6 = {°' 1 }- (2-32) 

Hence in a neighbourhood of z* = Z3, and z* = Z4 local solutions have the form [ 2.20) . 
Thus, we proceed in a similar way as in the proofs of previous lemmas. 
There exists only one e 6 £ = nf=i E; / such that d(e) G No- This is 

/llll 1 \ 
e= ( 2 - 33 ) 



4'4' 2 '2' 2 

and for it we have d(e) =0. Then we construct rational functions 



5 e- 

and check whether monic polynomial P = 1 satisfies differential equation l |2.24) for some 
values of parameters a and 5. This equation is satisfied only when equalities 

2a = 0, 1 + 2a 2 + a 2 S 2 = 0, a + <x 3 5 2 = 0, 

hold but this system of algebraic equations is not consistent. □ 

Lemma 2.6. If 5 = 1/a and 8 ^1, then the differential Galois grout) of equation [ 2.16) with 
r{z) given by ( |23l) is SL(2, C). 

Proof. In this case r given in 1 2.31) simplifies to the form 

3 3 1 l+3a -l+3« 

' ~ ~16(z-l) 2 ~ 16(z + l) 2 ~ 4(z-z 3 ) 2 ~ 16(a-l)(z-l) + 16(« + l)(z + l) 
a 3 a 



2(-l+a 2 )(z-z 3 ) 2(z-z 4 )' 

(2.34) 

where Z3 = a -1 and Z4 = — a -1 . Singularities Zi 2 = ±1 and Z3 are poles of the second 
order, and Z4 is a pole of the first order. The infinity Z5 = 00 is generically of degree 3, i.e., 
it is a regular singular point. Its degree changes only for excluded values of a, i.e., for 
a G {0, 1}. The difference of exponents at z* = Z3 is zero. Thus local solutions around 
this point have the form [ 2.20) , and the differential Galois group is either a subgroup of 
the triangular group, or it is SL(2, C). Thus we have to check only the first case of the 
Kovacic algorithm. 

The sets of exponents £, are now following 

El = E 2 = E3 = {^}' E4 = E 5 = {0,1}. 

Only for e = (1/4, 1/4, 1/2,0, 1) G nf=l' vame of d i e ) is an integer equal to d(e) = 0. 
Inserting rational function 



4 ft 



to 



E 



1=1 ' 



and P = 1 into differential equation [ 2.24) we obtain a(z + a) =0. This condition is ful- 
filled only for excluded value a = 0. Thus, the differential Galois group is not a subgroup 
of the triangular group, so it is SL (2, C ) . □ 
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Case when 3 = 1 corresponds to 7 = and now we prove the following lemma. 
Lemma 2.7. For 7 = differential Galois group of \2.16\ is SL(n, C). 

Proof. Using particular solution q> defined by | |2.30| one can easy check that the first and 
the second case of the Kovacic cannot occur but an analysis of the third case is very oner- 
ous. 

Fortunately in this case we can find another particular solution for which analysis of 
the differential Galois group of the corresponding variational equations is simpler. 

Namely Hamilton equation governed by Hamiltonian | |1.3) for m\ — 0, 7123 = m 2 and 
h — h have another particular solution lying on invariant manifold 

N= {(71/ 72, 7t lf K 2 ) I 72 =0 7T 1 =27r 2 }, 

We consider only the level c = 0. It is convenient to make the following non-canonical 
change of variables 



"71" 




~X\ 




72 


= A 


X 2 


A = 






V\ 




J 1 !. 




J/2. 





1 

10 

1 

1 1 

2 2 







(2.35) 



After this transformation equations of motion take the form 



x = JV*H, 



BIB , B 



H = H(Ax), 



(2.36) 



where x = \x\,x-i, y\,yi\ T ■ Hamilton function of the reduced system in these variables 
takes the form 

H = Wlyl + til{y\ + yl) - /i(cosx 2 (-4/ 2 y 2 (yi +yi) + h{yi -y 2 ) 2 cosx 2 ) 

+ 2(yi -y 2 )(2/ 2 y 2 + /i(yi + y 2 ) cosx 2 ) cos(xi +x 2 ) 

+ h (yi + y2) 2 cos(x! + x 2 ) 2 )] I l&{llm 2 {2 - cos(2x 2 ) - cos[2(x! + x 2 )])]. 



In these variables our system has invariant manifold 

N = {(x 1 ,x 2 ,yi,y 2 ) \x 1 =y 1 =0}. 
Equations of motion restricted to this manifold read 

yiiJl + ll +2/i/ 2 cosx 2 ) 



(2.37) 



x 2 



y2 



U\l\m2 



sin x 2 



y|(2(Zf + l\) cosx 2 + Z 1 Z 2 (3 + cos(2x 2 ))) 
4Z??2 m 2 sm3 x 2 



(2.38) 



These equations defines our particular solution and variational equations along it take 
the form 



(2.39) 



X{ 




"fl n 





«13 


" 




"Xi" 


x 2 




«21 


«22 


«23 


a 24 




x 2 






«31 





«33 









Y 2 




Ml 


fl 42 





a 44 _ 




y 2 
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where 

y z {l 2 + hcosx 2 ) 1 
Tlil^mi'smxi l 2 m 2 ) 



«31 



y 2 ({71 j + 8Z5) cosx 2 + k(10l 2 + 6Z 2 cos(2x 2 ) + h cos(3x 2 ))) 
16Z 2 Z 2 ?«2 sin 3 x 2 
_ y 2 (2 ( Zf + Z 2 ) cos x 2 + h l 2 (3 + cos (2x 2 ) ) ) 1 

" 2/2/2^2 sin 3 x 2 ~ fl23 "2lfm2' 

Z 2 + l\ + 2/^2 cos x 2 
2Zj Z|ffi2 sin 2 *2 

yl(10hl 2 cos x 2 + 4(Z 2 + 2Z 2 ) cos(2x 2 ) + Zi(6Z 2 cos(3x 2 ) + Zi (3 + cos(4x 2 )))) 

32Z 2 Z 2 OT2 sin 4 x 2 

y 2 {h + h cosx 2 ) 



a 2l 



a 2i 



2lil2tn 2 smx 2 

1 y\(23hl 2 cos x 2 + 4(Z 2 + Z 2 )(2 + cos(2x 2 )) + Z1Z2 cos(3x 2 )) 



«44 



2 16Z 2 Z 2 wz 2 sin 4 x 2 

y 2 (2 (Z 2 + Z 2 ) cos x 2 + hl 2 (3 + cos(2x 2 ))) 



2Z 2 Z 2 W22 sin 3 X2 

Thus we have the closed subsystem of normal variational equations 

Xj = flnXx + fl^Yx, Yj = flgiXx + fl33Vi- (2.40) 
This system can be rewritten as one second order linear equation 

To rationalise this equation we use transformation 

fi — >z = cosx 2 (f). (2.42) 
This transformation together with the following expressions on time derivatives 
d_ _ .d_ S_ _ 2 d 2 ..d_ 

dz ~ z d? dz 2 ~ {z) d? +z d? 

._ . . y 2 (Z 2 + Z 2 + 2ZiZ 2 cosx2) 2 y 2 (Z 2 + Z 2 + 2ZxZ 2 z) 2 

Z — — SIIII2I2 — ; / ( Z) 



2Z 2 Z 2 m 2 sinx 2 ' tt\l\m\(-\ + z 2 ) ' 



.2 . .. y^(Z 1 2 + Z| + 2Z 1 Z 2 z) 

Z = -COSX2X 2 -SinX2*2 = ~ ,,3,3 2/ 1 i 2T ' 

4ZjZ|wz|(-l +z 2 ) 

translates equation ( 2,41| into the following one 



X" + p(z)X' + ^(z)X = / X = Xj, (2.43) 

with rational coefficients 
z 1 

P ^F~2(z-z 3 )' 

& 1 a 2 a 2 2a 3 

+ 



1 z 2 4(z-z 3 ) 2 2(l+a) 2 (z-l) 2(a-l) 2 (z + l) (a 2 - l) 2 (z - z 3 ) ' 
Here z 3 = — (1 + a 2 ) / (2a) and a = Zi/Z 2 . Equation 1 2.43) has four singularities 



z 1/2 = ±1, z 3 = , z 4 = 00. 

2a 
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Using the transformation [ 2,15) to 1 243} we obtain its reduced form 

1 1 

w" = r{z)w, r(z) = -q(z) + -p'{z) + -p(z) 2 , (2.44) 

where the expansion of coefficient r is the following 

- - 7 - * 2 OL 1 _ 2o? 

r " " 16(z - z 3 ) " 2(a + l) 2 (z - 1) + 2(a - l) 2 (z + 1) " _ i)2( z _ Zs ) ^ ' 5 ' 

and order of infinity is 2 provided a ^0. 

At first let us assume that d ^ 1. To check the differential Galois group of equa- 
tion ( 2.44| we apply the Kovacic algorithm. In the first case sets £/ are following 

£ 1 = £ 2 = {0,1}, £ 3 = £ 4 = |^-i^4 +i x}- (Z46) 

Next, we select from the Cartesian product E — Ylf = i Ej those elements e = [e\, £2, £3, £4) 
for which 

3 

d(e)=e i -J^e t eN . (2.47) 
(=1 

We have two choices of e £ E, namely 

&=(o.o.\-4.\-4). & = (0,0,1+4. 1+4). (2.48, 

that give d(eW) = for z = 1,2. Then we construct rational functions 

and substitute into equation 1 2.24| setting also P = 1 . This equation is satisfied only for 
a = for both choices eW, z = 1, 2. 

In the second case of the Kovacic algorithm sets E, are the following 

£1 = £ 2 = {4}, E 3 = E 4 = |2-iV3,2,2 + i^} HZ = {2}. 

Thus we see that there is no element e = £3, £4) for that 

d(e) = ±(e4-pjj ef% (2.49) 
Similarly in the third case of this algorithm sets £, are following 

where n = 4, 6, 12. But we see that also there is no e = {e\,ei, £ 3 , £4) with e; £ E, given 
above for that 

Thus differential Galois group is whole SL(2, C). 

For a = 1 we have a confluence of singular points. In this case coefficient r simplifies 
to the following forms 

1 + 7z 

r ~ ~16(z-l)(l+z) 2 ' 
We note that in these cases equation < |2.44) becomes a Riemman P equation. With a help 
of the Kimura theorem 1 6 ] one can show that the differential Galois group of this equation 
isSL(2,C). 1 ' □ 
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The second part of our main Theorem 1 . 1 is an immediate consequence of Theorem 2 . 1 
and Lemmas |2 . 5\\2 . 6 1 and \2.7\ 



3 Conclusions 

The integrability analysis for the flail pendulum with parameters satisfying 

m 2 l 2 = m 3 h 



(3.1) 



was almost finished except the case when m\ = 0, m 2 ^ m 3 and l 2 ^= I3. The only 
integrable case was found for l\ — 0. 

There is an open question about integrability in the case when equality | |3.1) does not 
hold. Poincare sections suggest non-integrability for generic values of parameters. 



' ) SM 
• , 'Si ' 



; vS#? 









■'.''.".■■■■/./■■ ' .. 






, ■ 





*;;S , 



0.6 0.4 



0.2 OA 



72 



Figure 7: The Poincare sections for m\ = m 2 = m 3 = 1, Zi = Z3 = 1, l 2 — 2, 5T3 = c = 1 and 
£ = 0.042 (on the left) and m l = m 2 = 1, m 3 = 2, Z x = Z 2 = 1, Z3 = 3, n 3 = c = 1 and E = 0.00358 
(on the right). Cross-section plane 71 = 
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A Explicit forms of reduced equations and variational equa- 
tions 



The explicit form of vector field v = [v\, v 2 , v 3 , v±] in Section 2.1 defined by reduced 



equations of motion \2A\ after linear transformation \2.3\ is the following 

vi = 2l\l 3 m 2 m 3 (cl 3 + l 2 y 2 - h(yi +1/2)) - l\( m l + m 2 + m 3 ) {}\m 2 y 2 + hm 3 (l 3 y x ~ hyi)) 

+h{hhm 2 3 {l 3 y l - l 2 y 2 ) cos{x x - x 2 ) 2 + l 2 m 2 cos x 2 (l 3 m 3 (l 3 (c - 3yi -y 2 )+31 2 y 2 ) 
+hl 2 m 2 y 2 cosx 2 ) + l 2 m 2 m 3 cos(x 1 - x 2 )(l 2 (l 3 (c - y x -3y 2 ) + l 2 y 2 ) -h(-l 2 y 2 

+l 3 {yi + y 2 )) cos x 2 ))} / l\l 2 l 2 (—m 2 m 3 (mi +m 2 + m 3 ) + m 2 m 3 (m 3 cos(xx — x 2 ) 2 + m 2 cosx^)) 
v 2 = -l 2 (m 1 +m 2 + m 3 )(l 3 y 1 - l 2 y 2 ) + l\m 2 (cl 3 + \ 2 y 2 - l 3 {y x + y 2 )) + h(hm(hyi 
— hyi) cos(xx — x 2 ) 2 + l 2 m 2 (cl 3 + 2I 2 y 2 — l 3 (2y x +y 2 )) cosx 2 — l 2 m 2 y 2 cos(xi — x 2 ) 
x (l 2 + Zxcos^ 2 ))] / l 2 l 2 {-hm 2 (mi + m 2 + m 3 ) + l 3 m 2 (m 3 cos{x 1 - x 2 ) 2 + m 2 cosx 2 )) 

^3 = [{h{-hyi + h{-c + y\ +yi)) + hhy 2 cos{x 1 -x 2 ) + h{l 3 yi - l 2 y 2 ) cos x 2 ) 
x {h{l 2 {l 3 m 2 + l 2 (4m 1 +m 2 +4m 3 ))y 2 - 2(l 2 - l 3 )m 3 (-l 3 y 1 + l 2 y 2 ) cosxi) cosx 2 sinx;i 
-hih + h)(hm 2 y2 + 2m 3 (-l 3 y 1 + l 2 y 2 ) cosxi) cos(3x 2 ) sinxx + 2/|m 3 (-/ 2 y 2 + l 3 (-c + y x 
+y 2 )) cos(2x 2 )sin(2x!) + h{l 2 {l 3 m 2 - l 2 (4m\ +3m 2 + 4m 3 ))y 2 cosx 1 + (Z 3 yi - l 2 y 2 ) 
x (2/ 3 (2(m! + m 2 ) + m 3 ) + {-l 2 + l 3 )m 3 cos(2^))) sinx 2 + 2/ 2 (c/ 3 + l 2 y 2 - h(y 1 +y 2 )) 
x {-l 3 m 2 + l 2 m 3 cos(2x 1 )) sin(2x 2 ) + l x (l 2 + l 3 ){l 2 m 2 y 2 cosx 1 +m 3 (-l 3 y 1 + l 2 y 2 ) cos(2x 1 )) 

x sin(3x 2 ))] / \ 2 \ 2 \\{2m\ + m 2 + m 3 —m 3 cos[2(x 1 — x 2 )] — m 2 cos(2x 2 )) 2 

^4 = P(/ 2 (-/ 2 y 2 + / 3 (-c + yi +y 2 )) + l 1 l 2 y 2 cos(x 1 -x 2 ) +Z 1 (/ 3 y 1 - Z 2 y 2 )cosx 2 ) 
x(2wi 3 cos(x! - x 2 )(l 2 (-l 2 y 2 + l 3 (-c + y 1 +y 2 )) + l 1 (l 3 y 1 - l 2 y 2 ) cos x 2 ) +l 1 l 2 y 2 {2mi 

+m 2 + 2m 3 — m 2 cos(2x 2 ))) sin(x! — x 2 )] / l 2 l 2 l 3 (2mi +m 2 + m 3 — m 3 cos[2(x! — x 2 )] 
—m 2 cos(2x 2 )) 2 

(A.1) 
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Entries of matrix in variational equations ( |2,9) are the following 

c (2l 2 l 3 + h(h + h) cos x 2 ) sin x 2 

khih ~ h){-hm 2 - h{mi + m 2 ) + (l 2 + l 3 )m 2 cosx\)' 
2l 2 l 3 m 2 + l\{l 3 nii + (l 2 + l 3 )m 2 ) + l\l 2 {l 2 +3l 3 )m 2 cosx 2 



an = 

ai3 = 
a u = 
a n = 

«23 = 
fl 2 4 = 

«3i = 
«33 = 
«34 = 
«41 = 
«43 = 
#44 = 



l\l\m 2 {-l 2 m 2 - h(mi + m 2 ) + (l 2 + / 3 )m 2 cosx^) 
(l 2 -h)(21 2 h + h(h + h)cosx 2 ) 
^\hh{-hm 2 - l 3 {mi + m 2 ) + (l 2 + h)m 2 cosx^)' 

c '3('2 + h cos x l) sinx 2 
hh{h ~ h)(-hm 2 - l 3 (nii+m 2 ) + (l 2 + l 3 )m 2 cosx 2 )' 
-l\l 3 m 2 - l\{l 2 m 2 + h(mi + m 2 )) + lil 2 m 2 cosx 2 (-2l 3 + l x cosx 2 ) 
l\l 2 m 2 (-l 2 m 2 - Z 3 (m 1 +m 2 ) + (l 2 + l 3 )m 2 cosxl) 
1 (/3-?2)('2 + /icosx 2 ) 



l 2 l 3 m 2 l\l 2 {l 2 m 2 + l 3 {mi+m 2 ) - (l 2 + l 3 )m 2 cosx\)' (a.2) 
2c 2 {l 2 + l 3 )smx\ 
(h ~ h) 2 (-2l 3 mi - (l 2 + l 3 )m 2 + (l 2 + l 3 )m 2 cos(2x 2 ))' 

2c(l 2 + l 3 )(l 2 + h cos x 2 ) sinx 2 
~ hh{h - h){-2l 3 mi - (l 2 + l 3 )m 2 + (l 2 + / 3 )m 2 cos(2x 2 ))' 

2c(Z 2 + 1 3 ) sinx 2 
lil 3 (-2l 3 mi - {l 2 + l 3 )m 2 + (l 2 + l 3 )m 2 cos(2x 2 ))' 
2c 2 l 3 sin x 2 

(l 2 - l 3 ) 2 (-2l 3 mi - (l 2 + l 3 )m 2 + {l 2 + / 3 )m 2 cos(2x 2 ))' 

2cl 3 {l 2 + li cosx 2 ) sinx 2 
~ hhih ~ h){-2l 3 mi - {l 2 + l 3 )m 2 + (1 2 + J 3 )m 2 cos(2x 2 ))' 
2c sin x 2 

h(-2l 3 mi - {l 2 + l 3 )m 2 + (l 2 + l 3 )m 2 cos{2x 2 )Y 



B Linear second order differential equation with rational 
coefficients and Kovacic algorithm 

Let us consider a linear second order differential equation with rational coefficients 

w" + p(z)w' + q{z)w = 0, p(z),q(z) eC(z). (B.l) 

A point z = c E C is a singular point of this equation if it is a pole of p(z) or q(z). A 
singular point is a regular singular point if at this point p(z) = (z — c)p(z) and q(z) = 
(z — c) 2 q{z) are holomorphic. An exponent of equation | |B.l) at point z = c is a solution 
of the indicial equation 

p(p-l)+p p + qa = 0, p = p(c), q =q(c). 

After a change of the dependent variable z — » 1 /z equation dBTTb reads 



io" +P(z)o/ + Q(z)io = 0, 



We say that the point z = oo is a singular point for equation | |B.l) if z = is a singular 
point of equation | |B.2) . Equation ^B.l) is called Fuchsian if all its singular points (includ- 
ing infinity) are regular, see |T5ll4|. 

The Kovacic algorithm |7] allows to decide whether all solutions of equation I B.l I 
with rational coefficients p(z) and q(z) are Liouvillian. Roughly speaking, Liouvillian 
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functions are obtained from the rational functions by a finite sequence of admissible op- 
erations: solving algebraic equations, integration and taking exponents of integrals. For 
a formal definition, see e.g., (7). 

If one (non-zero) solution W\ of equation | |B.l) is Liouvillian, then all its solutions are 
Liouvillian. In fact, the second solution w^, linearly independent from W\, is given by 







/ A ex p 


[-!'} 


/ w\ 





w 2 = w \ 

Putting 

w = y exp 

into equation [ |B,1L we obtain its reduced form 

1 1 

y" = r(z)y, r(z) = -q(z) + -p'(z) + -p(zf 



(B.3) 



This change of variable does not affect the Liouvillian nature of the solutions. For equa- 
tion dB3), its differential Galois group Q is an algebraic subgroup of SL(2, C). The fol- 
lowing lemma describes all possible types of Q and relates these types to the forms of a 



solution of (B.3 1, see 171ITTI. 



Lemma B.l. Let Q be the differential Galois group of equation I B.3 1. Then one of four cases can 
occur. 

Case I Q is conjugate to a subgroup of the triangular group 

b ' 



T 



a 

a' 1 



a e C*,b e 



in this case equation ( |B.3| has an exponential solution of the form y = Pexp f to, where 
P 6 C[z] and to 6 C(z), 

Case II Q is conjugate to a subgroup of 



c 

c" 1 



ceC* u 



c 
-c" 1 



c e C* 



in this case equation ( |B.3| has a solution of the form y = exp J to, where to is algebraic 
over C(z) of degree 2, 

Case III Q is primitive and finite; in this case all solutions of equation < |B.3| are algebraic, thus 
y = exp f to, where co belongs to an algebraic extension of C(z) of degree n = 4,6 or 12. 

Case IV Q — SL(2, C) and equation < |B.3) to no Liouvillian solution. 

Kovacic in paper [7] have formulated the necessary conditions for the respective cases 
from Lemma [B.l I to hold. 

At first we introduce notation. We write r(z) E C(z) in the form 



r z 



s(z) 
f(z)' 



s(z), f(z) e C[z], 



where s(z) and f(z) are relatively prime polynomials and f (z) is monic. The roots of f (z) 
are poles of r(z). We denote X/ := {c £ C | t(c) = 0} and E := U {oo}. The order 
ord(c) of c 6 X/ is equal to the multiplicity of c as a root of £(z), the order of infinity is 
defined by 

ord(oo) :— deg t — deg s. 



Lemma B.2. The necessary conditions for the respective cases in Lemma B.l are the following. 
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Case I. Every pole ofr must have even order or else have order 1. The order ofr at oo must be 
even or else be greater than 2. 

Case II. r must have at least one pole that either has odd order greater than 2 or else has order 2. 

Case III. The order of a pole ofr cannot exceed 2 and the order ofr at oo must be at least 2. If the 
partial fraction expansion ofr is 

then A, = y/\ + 4a; £ Q,for each i, Y,j bj = and if 

! ' ;' 

then y/1 +4g £ Q. 

In IZ| Kovacic also formulated a procedure, called now the Kovacic algorithm, which 
allows to decide if an equation of the form |B.3l possesses a Liouvillian solution and 
to find it in a constructive way. Applying it, we also obtain information about the dif- 
ferential Galois group of the analysed equation. Beside the original formulation of this 
algorithm several its versions and improvements are known |3, 8 11]. 

Now we describe the Kovacic algorithm for the respective cases from Lemma 



B.l 



Case I 

Step I. Let T is the set of poles of r. For each c £ T U {oo} we define a rational function 
[\fr\ c and two complex numbers , cC^ as described below. 

(ci) If c 6 T and ord(c) = 1, then 

Wr\ c = 0, «+ = cc- = 1. 
(c 2 ) If c 6 r and ord(c) = 2, then 

Wf\c = o. 

If r has the expansion of the form | |B.4) , then 

ccf = X - (l ± (B.5) 

(C3) If c 6 T and ord c = 2v > 4 with k > 2 (only even orders are admissible in this case), 
then 



r 



(z-c) v (z-c) 2 ' 



(z-c) v+1 \(z-c) 1 
and 

= (±- + v 



2 \ a 

(<x>i) If the order of r at infinity is > 2, then 

[Vr]co=0, a+=0, = 1. 
(002) If the order of r at infinity is 2, then 

[Vr]oo = 0. 



19 



If the Laurent series expansion of r at oo takes the form 



then 



",_ A ) (B.6) 



(B.7) 



o± = l(±--v). (B.9) 



(003) If the order of r at 00 is — 2v < (necessarily even in this case), then 

Vr co =az v + ---+S, (B.8) 
is the indicated part of the Laurent series expansion of y/r at 00 and 

r-[^]i = bz v ^+o(z v - 2 ). 

Then 

2 V « 

Step II. For each family s = (s(c),s(oo)), c € T, where s(c) and s(oo) are either + 
either — , let we compute 

If d is a non-negative integer, then 

cv(z) = £ Lc) [V~ r ] c + + s (oo) [0^, (B.IO) 

is a candidate for <x>. If there are no such elements, equation \B.3\ does not have an expo- 
nential solution and the algorithm stops here. 

Step III. For each family from step II that gives d E No we search for a monic polynomial 
P = P(z) of degree d satisfying the following equation 

P" + 2cv(z)P' + (w'{z) +co{z) 2 - r(z))P = 0. (B.ll) 

If such polynomial exists, then equation (|B.3|> possesses an exponential solution of the 



form y = Pexp J 6, where 8 — co, if not, equation |B.3l does not have an exponential 
solution. 

Case II 

Let T be the set of poles of r. 
Step I. For each cer we define E c as follows 

(ci) If ord(c) = 1, then the set E c = {4} i 

(C2) If ord c = 2 and r has the expansion of the form < |B.4) , then 

E c := {2,2(l + Vl+4a), 2(l- Vl+4«)}nZ, (B.12) 

(c 3 ) If ord c = v > 2, then E c = {v} 

(001) If ord 00 > 2, then £ ro = {1,2,4} 

(00 2 ) If ord 00 = 2 and the he Laurent series expansion of r at 00 takes the form < |B.6) , 
then 

£co := {2,2(1 + VT+Iaj, 2(l- Vl+4a)}nZ. (B.13) 

(00 3 ) If ord 00 = v < 2, then £00 = {v} 
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Step II. We consider all families (e c ) ce ru{oo} with e c £ E c and at least one of the coordi- 
nates is odd. We compute 

and select those families (£c) c eTu{co} f° r wmcn d{e) is a non-negative integer. If there are 
no such elements, Case II cannot occur and the algorithm stops here. 
Step III. For each family giving d £ No we define 

to = to(z) = Je^7' ( b - 14 ) 

cer z L 

and we search for a monic polynomial P = P(z) of degree d satisfying the following 
equation 

P'" + 3loP" + (3a; 2 + 3co' - 4r)P' + (a?" + 307a;' + a; 3 - 4ra; - 2r')£ > = 0. (B.15) 



If such a polynomial exists, then equation | |B.3) possesses a solution of the form y 
exp J 6, where 



If we do not find such polynomial, then Case II in Lemma B.l 
Case III 

Step I. For each c £ T U {oo} we define a set E c of integers as follows 

(c\) If c is a pole of order 1, then E^_ = {12}. 

(C2) If ord c = 2 and r has the expansion of the form I B.4 \, then 



cannot occur. 



E c = vT + 4^|Jfc = 0,±l,±2...,±?l 

[ rw 2 J 



n z. (B.16) 



Here and below in this case m 6 {4, 6, 12}. 

(00) For c — 00 when the Laurent expansion of r at infinity is of the form < |B.6) we define 
independently on order 

£00 = <6+— VT + 4fl|Jfc = 0,±l,±2...,±^lnZ. (B.17) 
[ m 2 J 

Obviously it can appear that «co = 0. 
Step II. For e E E we calculate 

d(e) ■= H (coo - I . 



We select those elements e E E for which rf(e) is a non-negative integer. If there are no 
such elements, Case II cannot occur and the algorithm stops here. 

Step III. We consider all families (e c ) c eru{oo} with e c £ E c that give d(e) = n £ No . For 
each such family we define 

a; = co z) = — > . 

w I? ^ z — c 

Next we search for a monic polynomial P = P(z) of degree n satisfying a differential 
equation of degree m + 1 defined in the following way. Put P m = P. Then calculate P, 
for / = m, m — 1, . . . , 0, according to the following formula 

P/_! = -SP[ + [{m - i)S' - Sco]Pi -(m- + l)S 2 rP i+l , (B.18) 
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where 

s = U( x ~ c )- 



C6E' 



Then P_j = gives the desired equation for P. If such polynomial exists, then equa- 
tion |B.3 1 possesses a solution of the form y — exp J 6, where 6 is a solution of the equa- 
tion 



Qip. 

' 7 e' = o. 



If we do not find such polynomial, then we repeat these calculations for the next element 
m E {4, 6, 12}. If for all m such polynomial does not exist, then Case III in Lemma 
cannot occur. 



B.l 
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